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Cj ' We theoretically investigate collective phase synchronization between interacting groups of glob- 

,^rfH ' ally coupled noisy identical phase oscillators exhibiting macroscopic rhythms. Using the phase 

^J , reduction method, we derive coupled collective phase equations describing the macroscopic rhythms 

f^ ■ of the groups from microscopic Langevin phase equations of the individual oscillators via nonlinear 

Fokker-Planck equations. For sinusoidal microscopic coupling, we determine the type of the collec- 
tive phase coupling function, i.e., whether the groups exhibit in-phase or anti- phase synchronization. 
We show that the macroscopic rhythms can exhibit effective anti-phase synchronization even if the 
microscopic phase coupling between the groups is in-phase, and vice versa. Moreover, near the 
onset of collective oscillations, we analytically obtain the collective phase coupling function using 
center-manifold and phase reductions of the nonlinear Fokker-Planck equations. 
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(^ ' Systems of limit-cycle oscillators are used to model various rhythmic phenomena in natu- 

^^ , ral and artificial systems. When the interaction between the oscillators is weak, the system 

can generally be described as coupled phase oscillators. Qualitative understanding of syn- 
chronization phenomena, in particular the emergence of collective oscillations, has been 
achieved successfully through this type of reduced model. In the present paper, we de- 
velop collective phase reduction method for macroscopic oscillations arising from mutual 
^ . synchronization of many microscopic oscillators within a group, which is typical in the real 

world. We show that weakly interacting groups of oscillators, each consisting of globally cou- 
pled noisy identical phase oscillators and exhibiting collective oscillations, can be reduced 
to coupled equations for collective phases. This makes it possible to analyze the nature of 
synchronization between the collective oscillations in a closed way at the macroscopic level. 

I. INTRODUCTION 

Populations of coupled rhythmic elements can exhibit macroscopic oscillations through mutual synchronization [l|- 
|6|. The phase oscillator models have played important roles in theoretically analyzing their behavior, and the special 
class of models given by globally coupled phase oscillators in particular, was studied most intensively in the past 
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4l2j . Theoretical predictions based on such models have also been experimentally validated, e.g., in electrochemical 
oscillator systems |13l - [l7| and in discrete chemical oscillator populations [l8l - l20| . 

Recently, macroscopic synchronization between interacting groups of globally coupled phase oscillators exhibiting 
collective oscillations has attracted attention |2l| - |26{ . In most of the works so far, the macroscopic properties such 
as mutual entrainment between the groups have been analyzed through the microscopic individual phases. However, 
because we are interested in the macroscopic behavior of the collective rhythms exhibited by the oscillator groups, it 
should be much more convenient if each group of oscillators can be treated as a single macroscopic oscillator. Based on 
such consideration, we have developed collective phase reduction methods ^27-29], which provide us with the collective 
phase sensitivity of macroscopic rhythms of the oscillator group to weak perturbations. 

In this paper, we employ the notion of collective phase description |27H29J . and formulate a theory for weakly 
interacting groups of globally coupled noisy identical phase oscillators in a closed form at the macroscopic level. 
Specifically, we derive coupled collective phase equations from microscopic Langevin phase equations describing weakly 
interacting groups of globally coupled phase oscillators. A general formula that gives collective phase coupling functions 
from the microscopic phase coupling functions between the individual oscillators is obtained, and for the case with 
sinusoidal coupling, the types of the collective phase coupling function are determined as a function of the coupling 
parameters. Near the onset of collective oscillations, we can even analytically obtain the collective phase coupling 
function by the center-manifold and phase reductions. Based on the collective phase equations, we illustrate counter- 
intuitive phenomena in which two oscillator groups become anti-phase synchronized in spite of in-phase microscopic 
coupling between the groups, and vice versa (in-phase synchronization despite anti-phase microscopic coupling). 

In Ref. [30|, we considered a similar problem, namely, collective phase synchronization between two groups of 
globally coupled oscillators. The crucial difference is that we treat noisy identical phase oscillators in the present 
work, whereas we analyzed noiseless non- identical phase oscillators in Ref. [30|. Although these two situations look 
similar, they are essentially different physical systems (i.e., stochastic vs. deterministic) and mathematical treatments 
should be developed independently. Here we apply center-manifold reduction as well as phase reduction to nonlinear 
Fokker-Planck equations governin g th e oscillator groups, whereas we used the Ott-Antonsen ansatz in the analysis of 
the noiseless non- identical system \3()\. In both cases, despite the large difference in their mathematical structures, we 
obtain similar coupled collective phase equations describing macroscopic dynamics of the groups. Thus, the present 
paper and Ref. [30| are mutually complementary and together give deeper understanding of macroscopic collective 
phenomena. 

The organization of the present paper is the following. In Sec. [Hi we introduce a model of weakly interacting 
groups of globally coupled noisy phase oscillators and illustrate both effective anti-phase and in-phase collective phase 
synchronization between the groups by numerical simulations. In Sec. IIIIl we develop a theory that derives coupled 
collective phase equations from the microscopic model and determine the effective type of phase coupling between 
collective oscillations. In Sec. IIV[ we analytically obtain the collective phase coupling function near the onset of 
collective oscillations and discuss several important cases. In Sec. |Vl we discuss a relation to noise-induced turbulence 
in a system of nonlocally coupled oscillators. Concluding remarks will be given in the final section. 

II. THE MODEL AND ITS DYNAMICS 

A. Interacting groups of globally coupled phase oscillators 

We consider two interacting groups of globally coupled noisy identical phase oscillators described by the following 
model: 

N N 

fc=i fc=i 

for j — I,--- ,N and (cr, r) = (1,2) or (2,1), where 4>f {i) is the phase of the j-th oscillator in the cr-th group 
consisting of N oscillators and uj is the natural frequency common to all oscillators. The second term on the right- 
hand side represents the internal coupling between the oscillators within the same group, the third term represents the 
noise, and the last term gives the external coupling between the oscillators belonging to different groups. The internal 
phase coupling function r(0) is assumed to be in-phase, dT{(j))/d(f)\^=o < Q, namely, the oscillators within the same 
group tend to synchronize with each other. The external phase coupling function between the groups is described 
by Fo-T (</>)■ Characteristic intensity of the internal coupling within each group is scaled to unity, whereas that of the 
external coupling between the groups is given by e > 0. The noise Q"^ (t) is assumed to be white Gaussian |3ll433j . 



whose statistics are given by 

(^r W) = 0' (^r W^r^^)) = ^^skS.rSit - s). (2) 

The noise intensity is characterized by Z? > 0. When the external couphng is absent, i.e., e = 0, Eq. ([1]) has a 
critical noise intensity Dc below which phase coherent states are realized, namely, collective oscillations arise when 
< D < Dc [2]. In the following, we assume that e is sufficiently small and each group of oscillators exhibits stable 
collective oscillations. 

B. Phase synchronization between collective oscillations 

In the following numerical simulations, we assume that the phase coupling functions are sinusoidal (note however 
that our theory itself can be applied to general 27r-periodic phase coupling functions ^28.] ) . Without loss of generality, 
the natural frequency can be assumed to be zero, w = 0. The internal phase coupling function between the oscillators 
within the same group is given by 

r(0) = -sin(</. + a), l«l<|, (3) 

which is in-phase (attractive). In this case, the critical noise intensity is given by Dc = (cosq;)/2 as explained in 
Sec. IIVI The external phase coupling function is described by 

r,.(0) = -sin(0 + /3), |/3|<^, (4) 

which can be either in-phase (attractive) (|/3| < 7r/2) or anti-phase (repulsive) (|/3| > 7r/2) [2]. Introducing a complex 
order parameter A^^^t) with modulus R'''^\t) and phase 0('^^(t) through 

A(-) (t) = i?(-) (t)e'«''' ^'^ = - E ^'*^"' ^*^ ' (5) 

fe=i 

we can rewrite Eq. ([T]) with the sinusoidal coupling functions given in Eqs. ^ and (jlj as follows: 

^Y^ {t)=uj- i?(") sin (4"^ - ©'"' + a) + %/I? C]"^ (t) - ei?(") sin (0^"' - O^") + /?) . (6) 

Note that R^"^' quantifies the degree of synchronization and Q^°'' gives the collective phase of the u-th group. 

Focusing on weakly coupled collective oscillations, we carried out numerical simulations of Eq. ^ with Eq. ([5]) 
under the following conditions: the external coupling was assumed to be much weaker than the internal coupling, 
i.e., e — 0.01 ^ 1; we set D = Dc/2 = (cosa)/4 and a — Stt/S; the number of oscillators in each group was 
iV = 10000, which was sufficiently large to observe clear collective oscillations. We separately prepared two groups of 
phase oscillators exhibiting collective oscillations and used these states as the initial conditions of the simulations. 

In Fig. [TJa), evolution of the collective phase difference jG^^-* — 0*^^^! from almost in-phase synchronized state of 
the groups is shown. In spite of the in-phase external phase coupling condition between individual oscillator pairs, 
/3 = Stt/S, the collective phase difference jO^^^ — 0^^^| eventually approached tt, namely, the two groups became anti- 
phase synchronized after some time. Thus, Fig. [TJa) indicates that the collective phase coupling function between 
the group is anti-phase although microscopic external phase coupling functions are in-phase. In contrast. Fig. [Ijb) 
shows evolution of |0^^^ — 0'^'| from almost anti-phase synchronized state of the groups with anti-phase microscopic 
external phase coupling function, /3 — — Stt/S, which eventually became in-phase synchronized. 

Snapshots of the microscopic phase variables after the collective phase difference has reached the asymptotic value 
in Fig. [1] are displayed in Fig. [2] In Fig. [2](a), the two distributions of the oscillators are shifted by tt, indicating anti- 
phase synchronization between the groups. In contrast, the two distributions almost overlap in Fig. [2Ib), i.e., they 
are in-phase synchronized. Note that oscillators from different groups do not synchronize with each other, in other 
words, the collective phase synchronization between the groups is not due to complete synchronization of individual 
oscillators at the microscopic level. 

Thus, the type of the collective phase coupling functions can be effectively different from that of the microscopic 
external phase coupling functions, depending on the collective dynamics of the oscillators taking place in each group. 
We develop a theory that yields the collective phase coupling function from the microscopic model in the following 
sections. 



III. COLLECTIVE PHASE REDUCTION 

We derive coupled dynamical equations for the collective phase variables of the groups from the Langevin phase 
equations of individual oscillators through nonlinear Fokker-Planck equations and obtain a formula that relates the 
collective phase coupling function between the groups to the microscopic phase coupling function between individual 
oscillator pairs from different groups. Using them, we determine the type of the collective phase coupling function 
and explain the results of the numerical simulations in Sec. |IT1 

A. Nonlinear Fokker-Planck equations 

In the continuum limit, i.e., N — > oo, the Langevin phase equations ([!]) can be transformed into the following 
coupled nonlinear Fokker-Planck equations [2|, l2l|, [23, [2a] : 
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for {a,T) — (1,2) or (2,1). Here, f^°'^{(j),t) is the one-body probability density function of the individual oscillator 

phase 4> in the cr-th group, which is normalized as L d(j>f^'^' {<j>,t) = 1. The first two terms on the right-hand side 
represent internal dynamics of the cr-th group, and the third term represents weak interaction between cr-th group 
and r-th group. The complex order parameter of Eq. ([5]) is now expressed as 



A(-)(i)=i?(-)(i)e*«'"'W 
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When the external coupling between the groups is absent, each group of oscillators obeying Eq. ([7]) with e = exhibits 
collective rhythms under the condition < D < I?c 0, [U, [23, [23 ■ We assume that this situation persists even if e 
becomes slightly positive and the two groups interact with each other weakly. 

The collectively oscillating solution of the nonlinear Fokker-Planck equations Q without external coupling (e — 0) 
can be expressed as a steadily rotating wave packet on a periodic interval [0, 27r], 



/ 



i'^) 



/),t) = /o(¥>('^)), ^(-) = 0-e(-), e(-) = f]. 



(9) 



for cr = 1,2, where the /o(<y5) represents the steady functional shape of the wave packet, 9^'^-' is the location of the 
wave packet at time t, namely, the collective phase of the cr-th group, and fl is the collective frequency common to 
both groups (HHll- 



B. Collective phase equations 



Let us assume e = and focus on a single group. The group index a will be dropped for the moment. Inserting 
Eq. dH) into the nonlinear Fokker-Planck equation ([7]) with e = 0, we find that fo{>fi) (v? = 0— 6) satisfies the following 
equation: 
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Let u{ip^ t) represent small disturbance to the collectively oscillating solution and consider a slightly perturbed solution 
/{(pji) = fo{ip) + u{ip,t). Equation ([T]) with e = is linearized in u{ip,t), i.e., dtu{ip,t) — Lu(p,t), where the linear 
operator L is given by 



Lu (ip) = D-—rU {ip) + (O - w) -—u (ip) - -— .90 {^) u {ip) 
dip'' -'■- -^ - 
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Defining tlie inner product as 

u*{ip),u{ip) =/ dipu* {ip)u{ip), 

we introduce an adjoint operator L* of L by 

u* (ip) , Lu (if) = L*u* (ip) ,u{ip) 
The adjoint operator L* is explicitly given as 



(13) 



(14) 



tu* ip) = D-^u* ip) ~{n-co) ^u* ip) + go iv) ^u* (^) + j'^ V T ^ - p) /o (^') -^u* (^') . (15) 

In the calculation below, we need only zero eigenfunctions uo(</') of L and Uq{lp) of L* . Note that the right zero 
eigenfunction can be chosen as 



Luo [ip) = 0, uo ((/?) = -J- Jo if) , 

dip 



(16) 



which follows from differentiation of Eq. (jlOl) with respect to ip. The left zero eigenfunction is normalized as 

L*u*{v)^Q, [u*{^),uo{p)] = l. (17) 

Now let us introduce weak external coupling, i.e., we assume < e ^ 1 and treat the last term in Eq. ([7]) as 
perturbations. Using the phase reduction method [3, [23: [23] i "^^ can derive coupled collective phase equations from 
the nonlinear Fokker-Planck equations ([7]). Namely, we project the nonlinear Fokker-Planck equations ([7]) onto the 
unperturbed collectively oscillating solution as 

|(-e'^))^[.o(^-e(^)),|/(^)(,,.) 



-n-. 



(18) 



where we approximated f^'^'{4>, t) by the unperturbed solution foif ) and used that [uo{ip), /o(¥')] = ~^- Therefore, 
the collective phase equation takes the form 



e(^) = fi + er<xr (e^-^^-e^^) 



where the collective phase coupling Junction is given by 
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F^r (e(") - e(^)) ^ J\ipJ\ip' r,, [ip-ip' + e''^) - e(^') ko (ip) /o (^') 

for {a,T) = (1,2) or (2,1). The function ko{ip) is defined by 
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and normalized as 



ko i'f) = -Jo if) -j-u*o ((^) 
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dp ko (p) = / dpuo (p) Uo (<p) = 1, 
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(21) 
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which is the kernel function determining the collective phase sensitivity of the group as a convolution of the microscopic 
phase sensitivity |28| . 



C. The case with sinusoidal coupling 

When the microscopic external phase coupHng function T„r{4') is sinusoidal as given in Eq. (j4]), the collective phase 
coupling function also takes a sinusoidal form 

r^^(e) = -psin(e + <5) (23) 

because Eq. (POJ) is a double convolution of T^t {ip— ip' + 9) with feg {'P)fo if')- Here, the parameter p cos 5 determines 
the effective type of the collective phase coupling function; it is in-phase when pcosS > and anti-phase when 
pcosS < 0. This quantity can be evaluated from the following equation: 



. dFar (6) 

pcoso = - 



de 



dip I dip'T^ri-i^- •i^')koi(p)uoiip') ■ (24) 

e=o 



Now we examine the case D — Dc/2 — (cosq;)/4, which we considered in the numerical simulations shown in 
Fig. [H Typical functional shapes of /o(0), wo((/)), Uq(0), and fco(0) in this case are illustrated in Fig. [3l which were 
numerically obtained from the nonlinear Fokker-Planck equation. Details of the numerical method are described in 
Ref. [231 • From these functions, the dependence of pcoaS on a and (3 was numerically evaluated by Eq. (|24|) as shown 
in Fig. mja). The type of the collective phase coupling function is represented in Fig. IHJb), where the solid curves 
satisfying pcos6 = represent the borders between the in-phase and the anti-phase parameter regions. The two sets 
of parameter values used in Fig.[T]are also plotted in Fig.Hl^b). As can be seen, the set of parameters corresponding to 
Fig.[IJa) is in the anti-phase region, pcosS < 0, which yields effective anti-phase collective phase coupling between the 
groups. Similarly, the parameter set corresponding to Fig. [IJb) is in the in-phase region, pcosS > 0, yielding effective 
in-phase collective phase coupling. Thus, the collective phase reduction theory successfully explains the numerical 
results in Fig. [TJ 

IV. CENTER-MANIFOLD AND PHASE REDUCTIONS 

In this section, we analytically determine the collective phase coupling function at the onset of collective oscillations 
by applying phase reduction to amplitude equations obtained by the center-manifold reduction of the nonlinear Fokker- 
Planck equations. This method gives analytical results without recourse to numerical determination of the kernel and 
other functions for general microscopic phase coupling functions, though restricted to the vicinity of the onset of 
collective oscillations. 

A. Amplitude equations near the onset of collective oscillations 

We derive coupled amplitude equations that describe the macroscopic rhythms of the groups near the onset of 
collective oscillations. Expanding the 27r-periodic functions f^'^\(f),t), r(0), and Var{4') into Fourier series as 

f^^'^ (0, = ^ E fi^^ (^) ^'''*' //"' (*) = / ^"^ ^'"' ('^' ^) ^"'*' (2^) 

^— — oo 



r(0)= Y. r^^''*' ^' = ^J '^'^^ 



(0)e-''^ (26) 
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the coupled nonlinear Fokker-Planck equations ([7]) can be expressed as 






(28) 



When the noise intensity D is decreased below the critical value Dc in the absence of external coupling between 
the groups, e = 0, the uniform solution f''"\(j)^t) — l/(27r) of Eq. ([7]), corresponding to the incoherent state, becomes 



unstable. Equivalently, the trivial solution f^' — 1, f^' == (Z 7^ 0) of Eq. (f28| is destabilized and a pair of modes 
/j-; with critical non-zero wavenumbers ±Zc start to grow. From the linear part of Eq. ()28p . instability of the mode 
I occurs when Dl^ — ZlmF; < 0, namely, D < ImTi/l. Thus, the most unstable wavenumbers ±Zc are those that 
maximize ImTi/l. Generally, the fundamental harmonic components tend to be predominant in the phase coupling 
function, so that we obtain l^ — ±1 in most cases. 

Assuming Ic — ±1, we introduce a complex amplitude A^°''{t) of the fundamental harmonic modes of f^'^'{(j>,t) as 



/ 



2n 27r V / 



(29) 



where A'^'^^t) is the complex conjugate of A^°'^(i) — f_i{t). We are concerned with weakly coupled collective oscil- 
lations, and thus consider the external interaction as perturbations. Using the center-manifold reduction method [2], 
we can derive a pair of coupled complex amplitude equations from Eq. (|28p in the following form J21| : 



for {a, t) = (1, 2) or (2, 1), where the parameters are given by 

A* = -De — -D, -De = — Imr_i, rtc = ^ + Rcr_i + Fq + cFo-t.o, 

and 

-F_i(F_2+Fi) _ 

^ ~ 2ImF_i - zRcF_i + zF_2 ' ""-^ ' ' ' — i' 

See Refs. 0, lU, [23] for details of the derivation. We should note that Eq. ([5(1)) represents two coupled Stuart-Landau 
oscillators, each of which (i.e., A = {fi + if7c)^ — g\A\'^A) describes collective oscillations of the respective oscillator 
group. 



(30) 



(31) 



(32) 



B. Phase reduction of the amplitude equations 



Next, we derive coupled collective phase equations by reducing the coupled Stuart-Landau equations obtained above 
by assuming that the external interaction between the groups is sufficiently weak, i.e., e is small. When the two groups 
are uncoupled, e = 0, the limit-cycle solution Ao{0) of Eq. (j30]) is given by (a is dropped again for the moment) 
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The left and right Floquet eigenvectors of this limit-cycle solution associated with the zero eigenvalue can be written 
as 



u,ie) = '-^^^^ 



M .^.e^ 
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de \jReg ' '^^ ' y /I Reg 

where the inner product of Uq{Q) and U^{Q) satisfies the normalization condition 
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Though Eq. p4p is expressed in complex representation for the sake of convenience in analytical calculations performed 
below, they are equivalent to the known results [2|] . 

Now let us introduce weak external coupling as perturbations, i.e., we assume < e <C /-t ^ 1. Using the phase 
reduction method J2| , we can derive the collective phase equation ([T9|) from the amplitude equation (pO)) . Namely, we 
project the amplitude equation pop onto the unperturbed limit-cycle orbit as 



e^'") = Re 
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(36) 
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where we approximated A^"^^ by the unperturbed solution Aq{Q'^'^'>) and used that Ke[UQ{Q)Ao{Q)] = fl. Thus, the 
reduced equation is obtained in the form of Eq. (J19p . and the coUective phase coupling function is given by 



Fc 



(^QiT) _ e(r)^ ^ j^^ \^fj* (^Q{a)^j ^^^^^ (0^"^)] ■ (37) 



By inserting the expressions of Eqs. ([?!]) ■ (15^ . ([55|. and ([M|) into the formula Eq. ([57)) . the collective phase coupling 
function P^t-(0) can be analytically obtained, which takes a sinusoidal form 

r,, (8) = -p sin (e + 5) , pe**' = ^ . (38) 

Reg 

Note that we have not assumed that the external phase coupling function Tar{(l>) is sinusoidal so far. The sinusoidal 
collective phase coupling function arises because we assume that collective oscillations exhibited by the groups of 
oscillators are near the supercritical Hopf bifurcation point. 

When the phase coupling functions are given by the sinusoidal forms, Eqs. ^ and (jlj, the parameters of Eqs. pip , 
and (1551) can be calculated as 

cos a „ sin a „ 3 sin a D tan a 

D,^^-, n, = uj-^-, n = oj +^-^^ (39) 



and 



^ d,, = \e-^P. (40) 



4 cos a — 2isina 2 

Inserting Eq. ^^ into Eq. ([35]), we obtain 

1 i 

pe"-^ = - (2 cos ;3 - tan a sin /3) + - (2 sin /3 + tan a cos 0) . (41) 

Therefore, the type of the collective phase coupling function is analytically found from the following quantity: 

p cos (5 = - (2 cos /3 — tan a sin /3) . (42) 

Reflecting the symmetry of the original model of Eq. ^ with respect to (a,/3) -^ — (a,/3) and </> -^ —0, Eq. ((31]) 
is symmetric about the origin in the a-P plane. The type of the collective phase coupling function at the onset of 
collective oscillations, i.e., D — D^, is represented in Fig.jS] which is very similar to Fig. |3J 

C. Several important cases 

Here, we consider three special and important cases of the collective phase coupling functions derived for the 
sinusoidal microscopic phase coupling functions, Eqs. ([3]) and (|3]), at the onset of collective oscillations. 

(i) The first case is a = 0, which indicates that the internal phase coupling function within the same group is 
antisymmetric. Inserting a — into Eq. (jH]), we obtain the following result: 

pe'' = ^e^^ (43) 

so that pcosS — (cos/3)/2. Thus, the collective phase coupling function has the same type as the microscopic external 
phase coupling function. The internal phase coupling function does not affect the type of the collective phase coupling. 
A similar scenario has been encountered in different models 28-30]. 

(ii) Several special values of the microscopic external coupling phase shift (3 comprise the second case. Inserting 
/3 — 0, ±7r, ±7r/2 into Eq. (jUJ, we obtain the following results: 

/3 = 0, pe*^ = i + ltana, (44) 

/3 = ±^, pe'' = -^-jtana, (45) 
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pe = =F-tana± -. 



(46) 



For antisymmetric external interactions, i.e., /3 — 0,±7r, the type of the collective phase coupling function coincides 
with the microscopic external coupling and is not affected by the type of the microscopic internal coupling phase shift 
a. In contrast, for symmetric external interactions, i.e., /3 = ±7r/2, the type of the collective phase coupling function 
is solely determined by the internal coupling parameter a, which can be either in-phase or anti-phase. 

(iii) The third case is /3 = a, namely, when the external coupling has the same phase shift as the internal one. 
Inserting (3 = a into Eq. (I4ip . we obtain the following result: 



■ r 1, . > 3sma 

pe = — 2 cosa — tana sma + i • 

4^ '4 



(47) 



Note that |/3| — \a\ < tt/2 in this case, namely, both internal and external coupling functions are in-phase. The type 
of the collective phase coupling function is anti-phase when tan^ a > 2. As we discuss below, this condition is the 
same as the condition for noise- induced turbulence in nonlocally coupled phase oscillators [27j . 



ON NOISE-INDUCED TURBULENCE 



Finally, we briefly discuss the relation between "effective anti-phase coupling" and "noise-induced turbulence" . In 
this section, our arguments do not assume that collective oscillations are near the onset. In Ref. [231, '^^ considered 
a system of nonlocally coupled noisy phase oscillators described by the following model: 



d_ 



(j>{r,t)=uj+ / dr' G {r - r')r {(t) {r,t) - (j> {r' ,t)) + VD ^ {r,t) , 



(48) 



where (j){r, t) represents the phase field of spatially extended oscillatory media, G{r) is a nonlocal kernel function that 
decays with the distance |r|, r(0) is the phase coupling function, S.{r,t) represents spatiotemporally white Gaussian 
noise, and D is the noise intensity. 

The Langevin phase equation (|48p can be transformed into a nonlinear Fokker-Planck equation in the following 
form: 



d 



2tt 



d0' r (0 -0') /(</-', r-,t) I /(0,r-,t) 

d0'r(0-</>'){G2V2/(0',r-,t)}/(0,r,i) 



92 



Here, we have expanded the nonlocal coupling term as 

/oo 
dr' G{r~ r') f {4>' , r' , t) = ^ GsnV^"/ (0', r, t) , 
^ — n 



(49) 



(50) 



n=0 



where 6*2™ is the 2n-th moment of G(r). 

The space-dependent order parameter A{r,t) is defined by 



A (r, t) = R (r, t) e'®^""'*) = f dr' G {r - r') f 



■e'*/(0',r',i). 



(51) 



where Q{r,t) can be considered as the space-dependent collective phase, f{(j),r,t) — fa{(t> — 'd{r^t)). Applying the 
phase reduction method to Eq. P^ . we obtained the following collective phase equation 



d_ 
di 



e{r,t)^n + j^v^e (r, t) + /2(ve (r, t)y + 



where v and /x are coefficients. In particular, the collective phase diffusion coefficient v was given by 



^27r <.27r 

i' = -G2 / dip dip' r ((^ - ip') ko (ip) uo {(fi') , 

Jo Jo 



(52) 



(53) 
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which can be negative and then induce spatiotemporal chaos (turbulence) . Details of the definitions and the derivations 
are given in Ref. 27]. 

Now, it is clear that Eqs. ([7|), ()19|) . and (|24|) describing interacting groups of globally coupled noisy phase oscillators 
are similar to Eqs. (H^ . ([5^ . and ([5^ describing a system of nonlocally coupled noisy phase oscillators. When the 
external phase coupling function is the same as the internal one, i.e., Tcrif) — r('p), Eq. (j24p is essentially equivalent 
to Eq. ([53|) . Namely, the instability condition, — F'(0) < 0, for in-phase collective synchronization between two groups, 
which gives the anti-phase condition for the sinusoidal coupling, coincides with the instability condition, P < 0, for 
spatially uniform solutions of the collective phase equation. Therefore, the phase diagram plotted in Fig. [5] using D 
and a with /? = a is the same as that we obtained for noise-induced turbulence in Ref. [27| . 

The above situation for collective oscillations at the macroscopic level is in parallel with the classical problem 
for phase oscillators at the microscopic level, in which the instability condition for in-phase synchronization of two 
coupled phase oscillators coincides with the instability condition for spatially uniform solutions of the phase diffusion 
equation [2]. 

VI. CONCLUDING REMARKS 

In the present paper, we considered two weakly interacting groups of globally coupled noisy identical phase oscillators 
undergoing collective oscillations. To analyze them, we adopted the idea of collective phase description, namely, we 
treated the collective oscillations of each group as a single macroscopic phase oscillator. We developed a theory 
that derives the collective phase coupling function between the groups from the microscopic external phase coupling 
function between individual oscillator pairs belonging to the different groups. Based on this theory, we illustrated 
counter-intuitive situations in which the two groups become anti-phase synchronized despite in-phase microscopic 
coupling, and vice versa. We also developed a theory that gives explicit analytical expressions of the collective phase 
coupling functions near the onset of collective oscillations. A complete phase diagram in the case of the sinusoidal 
internal and external coupling functions is summarized in Appendix. 

In our companion work [30| , we considered two weakly interacting groups of globally coupled noiseless non-identical 
phase oscillators and discussed their collective synchronization properties. The strong similarity in results between 
the two types of systems, one stochastic and the other deterministically random, is remarkable, while the theoretical 
methods employed are completely different between them. In particular, we found the same counter-intuitive phe- 
nomena, namely, the disagreement of the types between the collective phase coupling function and the microscopic 
external phase coupling function. 

The notion of collective phase description is convenient and powerful in analyzing complex macroscopic rhythms 
arising from systems of interacting microscopic dynamical elements. Further development of the theories will pro- 
vide useful viewpoints to understand various complex rhythms in real-world systems, in particular, their functional 
meaning. 

Appendix: phase diagram for sinusoidal coupling 

We here present a complete phase diagram for the case of the sinusoidal internal and external coupling functions. 
The type of the collective phase coupling function is found from pcosS given in Eq. ([M]) . which was numerically 
evaluated for /3 G [— 7r,7r] in the parameter region of D/Dc E [0.1,1.0) and a G [— a7r/2, a7r/2] with a = 0.9. In 
addition, we used the analytical formula equation (H^ on the Hopf bifurcation line, i.e., D = Dc = (cosa)/2. Phase 
diagrams in D and a with several values of /? are displayed in Fig. [71 
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FIG. 1: (Color online) Time evolution of collective phase difference |B'^' — 0'^'|. (a) Effective anti-phase collective synchroniza- 
tion with microscopic in-phase external coupling, /3 = Stt/S. (b) Effective in-phase collective synchronization with microscopic 
anti-phase external coupling, /3 = — 57r/8. The other parameters are e — 0.01, uj = 0, D — Dc/2 — (cosQf)/4, and a = Stt/S. In 
numerical simulations, the number of oscillator in each group is A'' = 10000. 
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FIG. 2: (Color online) Snapshots of the asymptotic states of the individual oscillators in Fig.[T] Only one in every five oscillators 
is plotted. Open circle (o) and plus (+) respectively indicate oscillator of group (1) and that of group (2). (a) Effective anti- 
phase coupling with microscopic in-phase coupling, /3 = Stt/S. (b) Effective in-phase coupling with microscopic anti-phase 
coupling, /? — — Stt/S. 
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FIG. 3; (a) Distribution function /o(0). (b) Right zero eigenfunction uo(4>)- (c) Left zero eigenfunction Uq{(J)). (d) Kernel 
function ko{<j))- Parameters are D — Dc/2 = (cosq)/4 and a = Svr/S, where order parameter amphtude is _R = 0.65328. 
/3 = 3tt/8 gives pcosS = —0.32264, whereas l3 = — 57r/8 gives pcos^ = 0.32264. 
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FIG. 4; (Color online) Effective type of phase coupling between collective oscillations with a G (— 7r/2,7r/2), /3 G [— 7r,7r], and 
D = Dc/2 = (cosa)/4, which is numerically evaluated by Eq. ([24}. (a) Dependence of pcos5 on a and /?. (b) The solid curves 
are determined by pcos(5 = 0. The filled circle (•) indicates a = /3 = Stt/S corresponding to Figs.[Ha) and[2ja). The cross (x) 
indicates a = Stt/S and /3 = — 57r/8 corresponding to Figs. [Hb) and[2jb). 
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FIG. 5: (Color online) Effective type of phase coupling between collective oscillations in a and /3 with D — Dc = (cosq)/2, 
which is analytically given by Eq. (|42p . i.e., pcosS — (2cos/3 — tanQsin/3)/4. (a) Dependence of pcosS on a and /?. (b) The 
solid curves are determined by p cos 5 — 0. 
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FIG. 6: (Color online) Effective type of phase coupling between collective oscillations in D and a with /3 = a. (a) Dependence 
of pcosS on D and a. (b) The solid curves are determined by pcosS = 0. 
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FIG. 7: (Color online) Phase diagram in D and a with (a) (3 = 0, (h) (3 = Stt/S, (c) /3 = 7r/2, (d) /3 = Svr/S, (e) ^ = tt, (f) 
^ = -57r/8, (g) /3 = -7r/2, and (h) /3 = -Stt/S. 



